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Abstract 

The oblique collision between an elastic disk and an elastic wall is numerically studied. We 
investigate the dependency of the tangential coefficient of restitution on the incident angle of 
impact. From the results of simulation, our model reproduces experimental results and can be 
explained by a phenomenological theory of the oblique impact. 
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I. INTRODUCTION 



Collisions are common phenomena in nature. In macroscopic scales, we often see colli- 
sions of balls in sports such as the baseball and the billiard. In such collisions, the initial 
kinetic energy of material dissipates into internal degrees of freedom. A part of total en- 
ergy is distributed into the translational motion and the rotational motion, while the other 
part is dissipated as elastic vibration, sound emission, heat, etc. As a result, collisions of 
macroscopic material are always inelastic. 

Inelastic collisions play important roles in granular materials^, Characteristic be- 
haviors of granular material come from inelastic collisions among particles. The Distinct 
Element Method (DEM) is a standard method of simulation for the granular materials^. 
DEM contains some phenomenological parameters such as the Coulomb's coefficient of fric- 
tion, dashpots, and so on. Nobody can determine such the viscoelastic parameter from 
the first principle. Even the determination of the simplest parameter, the coefficient of 
restitution (COR), is difficult. 

COR is a familiar parameter which is introduced in text books of the elementary physics. 
The normal COR is defined by the ratio of the normal components of the collision velocity 
before and after collision. Figure ^ is the schematic figure that a sphere is colliding with a 
stationary wall with initial velocity of its center of mass, v. The prime denotes post-colliding 
quantities. The coefficient of normal restitution e is defined by 

v c • n = -ev c • n, (1) 

where v c and are respectively the velocity of the contact point before and after the 



collision, e is assumed to be < e < 1. Historically, COR was first introduced by Newton 
Though many text books of elementary physics state that e is a material constant, man 
experiments and simulations show that e decreases as the impact velocity increases ja, 0, 
0, hJ llll ? ] except for the idealized situation of a collision between two rods, in which e 



may be determined by the ratio of length of the colliding rods|l2l. ll3L [la, [l6(- Recently, 
Louge and Adams have reported that e can exceeds unity in an oblique impact The 
mechanism of the exotic behavior has not been clarified yet. Thus, study of collisions is not 
matured and has room for development though the subject itself is familiar even for high 
school students. 
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In addition to e, the tangential COR (3 is also important to characterize oblique impacts, 
where (3 is defined as 

v ; • t = -pv c • t, (2) 

where and t are the post-collisional velocity at the contact point after collision and the 
unit tangential vector, respectively. (3 is a function of the angle of incidence 7 which is 
defined as 7 = arctan(t>t/t> n ) with v n = v r • n and v t — v c • t and is believed that possible 



values of (3 lie between -1 and 1 



18, 



The aim of our study is to investigate (3 in detail from numerical simulation. We study 
the relation between (3 and the angle of incidence in oblique collision in this paper. The 
organization of this paper is as follows. In §2, we will review the current status of the 
researches about inelastic collisions. In §3, we introduce our numerical model and setup of 
the simulation. Section 4 is the main part of this paper where we summarize the results of 
our simulation and compare them with the theoretical outcome. Section 5 is devoted to the 
discussion of our results. The last section is the conclusion remarks. In Appendix [21 we 
will give a description of the theory of the oblique impact. In Appendix El the derivation 
of Poisson's ratio of the square lattice system will be introduced. 

II. REVIEW 

In this section, we review current status of the study of inelastic collision. At the end of 
the 19th century, it became clear experimentally that COR depends on the relative velocity 

. In the normal collision of spheres, it is believed that COR 
obeys 1 — e a v when impact velocity v is low enough while e oc u _1//4 when v exceeds 
the critical value of plastic deformations. In fact, Sondergaard et al^ performed an impact 
experiment of ball bearings and glass spheres on lucite or aluminum plates and confirmed the 
dependency on the impact velocity as e oc v in the high speed impact. Johnson explained 
this dependency using the dimensional analysis His analysis was consistent with the 
experimental data. Kuwabara and Kono'll] investigated COR for low speed impact and 
derived the theoretical expression 1 — e a v 1 ^ which is consistent with their experimental 



results. Recently, some papers a 



with Kuwabara and Kono 



Dout quasi-static theory are published which are consistent 



2j,i2al2 



22]- 



On the other hand, Gerl and Zippelius performed simulation of a two-dimensional collision 
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of an elastic disk with an elastic wall [281] . The present authors performed two-dimensional 



simulations and confirmed that elastic models including 
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;hat of Gerl and Zippelius is not 



3l|. 



appropriate to characterize the quasi-static region [29, 

To characterize the oblique collision, Walton introduced three parameters: the coefficient 
of normal restitution e, the coefficient of Coulomb's friction fi Q , and the maximum value 
of the coefficient of tangential restitution /3o|l8|- Experiments have supported that his 
characterization adequately ca ptu re the essence of binary collision of spheres or collision of 



a sphere on a flat plate 



P 







20, 21 



221. Walton derives 



-l + /i (l+e) cot7 1 + 



mR 2 



00 



(7 > 7o) 
(7 < 7o), 



(3) 



where 70 is the critical angle, and m, R, and I are mass, radius and moment of inertia 



of spheres respectively. Experimental results are consistent with eq. 
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34, 



19, 20, 21, 2 



and estab 



ished 



Meanwhile, Maw et al. extended the Hertz theory of impact [32 

the theory of the oblique impact to be consistent with their experimental results |3fi 
The theory by Maw et al. has several advantages: (i) The physical mechanism is included 
in their theory but not in Walton's argument, (ii) Anomalous behavior near 7 = can be 
described in their theory. The disadvantage of their argument is that we cannot summarize 
the result of theory in a concise form as in eq.(|3*)l. Thus, we still use both Walton's expression 
and the theory by Maw et al. The details of the theory by Maw et al. and its application 
to our results is summarized in Appendix [X] 



III. OUR MODELS 



In this section, let us introduce our numerical model. Our numerical model consists 
of an elastic disk and an elastic wall (Fig. |2J). Both of them are composed of randomly 
distributed 1600 mass points. We use the Delaunay triangulation algorithm to connect all 



mass points with nonlinear springs 
points is described as 



381 ] . The spring interaction between connected mass 



V(x) 



— k a x + —kbX , 



(4) 



where x is a stretch from the natural length of spring, and k a and kb are the spring constants. 
In most of simulations, we adopt k a = 1.0 x mc 2 /R 2 and kb = 1.0 x 10~ 3 mc 2 / R A , respectively. 
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The width of the wall is 4 times as long as the diameter of the disk while the height of the 
wall is same as the diameter of the disk. Two sides of the wall are fixed. 

The interaction between the disk and the wall during a collision is introduced as follows. 
Figure El is the schematic figure of the interaction of surface mass points of the disk and the 
wall. When the distance / between the lower edge of the disk and the surface of the wall 
becomes less than the cutoff length which is the mean value of natural lengths of all springs, 
the surface particles of the disk feel the repulsive force, F(l) = aVo exp(— al)n, where a is 
300/i?, Vo is amc 2 R/2, m is the mass of the particle, R is the radius of the disk, c = yVE/p, 
E is Young's modulus, and p is the density, n is the normal unit vector to the surface. The 
reaction forces applied to the two points of the surface of the wall (point 1 and 2) are decided 
by the balance of the torques as F 1 (l) = -F(Z)n/(l + /i/Z 2 ) and F 2 (7) = -F(l)n/(l + l 2 /h), 
where k(i = 1, 2) is the distance between the point p and the point i (see Fig. |3J). 

In this model, roughness of the surfaces is important to make the disk rotate after col- 
lisions. To make roughness, at first, we generate normal random numbers whose average 
value is and then make the initial position of particles on surfaces of both the disk and the 
wall deviate with them. We choose the standard deviation of the normal random numbers 
5 as 5 = 3 x 10~ 2 i?, where R is the radius of the disk. All of the data presented here are 
obtained from the average of 100 samples in random numbers. 

For random lattice model, it is impossible to determine Poisson's ratio v and Young's 
modulus E theoretically. When we determine Poisson's ratio v and Young's modulus E of 
this model, we add the viscous force term in equation of motion which is proportional to 
the relative velocity of two connected mass points. By stretching the band of random lattice 
and calculating the ratio of the strains in the vertical and horizontal directions to the force 
and the ratio of the strain to the force when the vibration stops, we can obtain Poisson's 
ratio and Young's modulus. Figure Ufa) and Efb) are snap shots of the bands of random 
lattice made of 348 mass points before and after adding the force F = 3.0 x 10 2 mc 2 /R, 
respectively. We change the force from 2.0 x 10 2 mc 2 /R to 3.0 x 10 2 mc 2 /R and average 10 
samples of results to obtain v = (7.50 ± 0.11) x 10" 2 and E = (9.54 ± 0.231) x 10 3 mc 2 /R 2 , 
respectively. 

For comparison, we introduce other two lattice models for elastic disks: triangular lattice 
and square lattice disk(FigEJ). To investigate the effect of the structure of the disk, the wall 
is same as random lattice model. In order to remove anisotropies of lattice structure, we 
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put an exterior layer which is same as random lattice disk around the disk and bind all of 
them using the Delaunay triangulation algorithm. The triangular lattice disk is made by 
replacing the internal structure of the random lattice disk with the triangular lattice. Total 
number of mass points is same as that of random disk. Poisson's ratio and Young's modulus 
of the triangular 
limit respectively 

the random lattice disk with the square lattice and connecting all the mass points by the 
Delaunay triangulation algorithm. We introduce two spring constants: k a = k% for nearest 
neighbor interaction and k„. = k 2 for next-nearest neighbor interaction. In the continuum 
limit, Young's modulus E\i$ and Poisson's ratio v (see Appendix^ of the square lattice 



attice can be calculated theoretically as 1/3 and 2k a / v3 in the continuum 
3t|. The square lattice disk is made by replacing the internal structure of 



are expressed as 



1 ki + k 2 k\ — 2k 2 2 2 



E k 1 (k 1 + 2k 2 ) hk 2 



+ L u n * n y> ( 5 ) 



k\ + {k\ - Akl)n 2 x n 2 



(6) 



k^ + k^ + tkl-AkDnlnf 
where n x and n y are the unit normal vectors horizontal and vertical to the collisional plane. 
The derivation of eqs.jJHJ) and © is presented in Appendix [HI To recover the orientational 
symmetry and introduce roughness on the surface of disks, we introduce one-layer random 
lattice on the surface of the disks. 

We scale the equation of motion for each particle using the radius of the disk R as the 



scale of length and the velocity of elastic wave c = E / p as the scaling unit of velocity. 
As the numerical scheme of the integration, we use the fourth order symplectic numerical 
method with the time step At ~ lO~ 3 R/c. 



IV. RESULTS 



In this section, we explain the results of our simulation. The angle of incidence 7 is ranged 
from 5.7° to 80.5° while the normal component of velocity is fixed as 0.1c. The colliding disk 
has no internal vibration and rotation at release time. In order to eliminate the effect of the 
initial configuration of mass points, we prepare 100 samples of disk as the initial condition 
by using 100 sets of normal random numbers and average data of all samples. 

Figure |H1 shows the relation between cot 7 and the coefficient of tangential restitution (3. 
In this figure, cross points are the result of the impact between random lattice disk and wall, 

6 



and broken lines are eq.(|3*J). In eq.(J3J) we use the value e = 0.8 which is the approximate 
mean value of e in the range 2.5 < cot 7 < 6 in Figd The result of simulation shows that 
fto is 0.56 and /io is 0.18 which are close to the values observed in experiments of three 
dimensional impacts jigL IscJ. Thus, we re pro duce experimental tendencies of the oblique 



collision with the random lattice model 



20j | . Stars are the results of random lattice disk 



without roughness on the surface, in which (3 is close to —1. From this result, one can see 
that roughness on the surface is important for the rotation of the disk after collision. Plus 
points in FigJS] are the result of the triangular lattice model where the orientation of initial 
disks is same as that in FigED^a). In this model, j3 takes negative values in all range of the 
angle of incidence. This means that the disk made of triangular lattice is easy to slip on the 
surface. In addition, the result strongly depends on the initial orientation of the disk. Thus, 
the model of triangular lattice is inadequate to reproduce the tendency of experimental data. 

Figure d shows the relation between cot 7 and e. Although e is expected to be a constant 
because the normal velocity of the disk is fixed, COR depends on 7. In particular, for small 
cot 7, e decreases as cot 7 decreases. We will discuss this behavior in the later discussion. 

nnn 

j which was consistent 

with experimental data|20j, |2l|, 



22J, |36j, |37[ • According to their theory, all the region of the 



angle of incidence can be divided into three regimes. For each regime, (3 can be expressed 
as 



1. 1/fir] 2 < cot 7: 



(3 = — coscjix — /i^-e 

Pz 



l + cos( — + -(l-e- 1 ) 

& — 



cot 7, 



(7) 



2. Px/PzK 1 + e) < cot 7 < 1/fir] 



2. 



(3 = — COSLu(t 3 — t 2 ) — fi—[cosLu(t 3 — t 2 ) — cos Vtt 2 coso;(t3 — t 2 ) 

Pz 

H sin Vtt 2 sin ^(£3 — t 2 ) + e + cos Vtt^} cot 7, 



(8) 



3. cot 7 < P x /(3 z fi(l + e): 



(3 = -1 + n^- (1 + e) cot 7, 

Pz 



(9) 



where fi is the coefficient of friction, 77 is the constant dependent on Poisson's ratio defined 
in eq. ()A5|) . (3 X and (3 Z are constants calculated from mass, radius, and radii of gyration of 



material as (3 X = 3.02 and /3 X = 1 as shown in eq. (jA2|) . Q and u are respectively ir/2t c 
and (7r/2r]t c )^(3 x / (3 ZJ where £ c is a duration of a collision, ti determined by eq. (jA10|) is the 
transition time from initial stick motion to slip motion. t$ determined by eqs. (|A21j) and 



(|A20|) is the transition time from slip motion to stick motion. By calculating (3 at each value 
of cot 7 and interpolating them with cubic spline interpolation method, we can draw the 
theoretical curve. 

We compare the result of simulation of the oblique impact using the random lattice model 
with the theoretical curve(Fig|HJ). Here we used 77 = 1.015, e = 0.8 which is an average value 
of COR in FigEl and \x = 0.18 which is decided by comparing the slope in the small cot 7 
region with eq.(|3|). It is found that the result of random lattice model is consistent with 
the theory especially in small cot 7 region. In the intermediate region, the agreement of the 
data with the theory is worse than that of other regions. This tendency can be seen in some 



experimental results 



20 



22 



37]. 



Theoretical result by Maw et al. suggests that Poisson's ratio is not a crucial parameter 
to determine /3, while Poisson's ratio of the triangular lattice is much larger than that of the 
random lattice. To confirm that anomalous behavior of the triangular lattice comes from 
the specific lattice structure we simulate the collision by using the square lattice model. By 
changing the value of spring constants of square lattice disk and controlling Poisson's ratio, 
we investigate the dependency of (3q on Poisson's ratio. /3q are 0.49 and 0.51 when v = 0.1 
and v = 0.3, respectively. From these results, we confirm that Poisson's ratio is not a crucial 
parameter for (3q. 



V. DISCUSSIONS 



Here, we discuss the results of our simulation. We change the number of mass points of 
random lattice model and investigate the dependency on the system size. As the number of 
mass points becomes larger, there is a tendency for a graph to be flattened in the region of 
large cot 7. It can be seen as follows. When the model is composed of many mass points, 
irregularity of the surface of the random lattice diminishes as the size of the disk increases. 
As a result, (3q can take the stable value in the large cot 7. When the number of mass points 
is larger than 1600, all the mass points in the rectangle cannot be connected by the Delaunay 
triangulation algorithm jsSj]. Hence, our results are restricted to the case with the 1600 mass 
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points as the maximum value. 

We also investigate the influence of roughness of the surface. When the standard deviation 
5 takes 2.15 x 10 _2 i?, (3 increases monotonously as cot 7 increases. When S takes 3.0 x 10 _3 i?, 
(3 approaches the stable maximum value (3q = 0.56. For larger 5, the surface of the rectangle 
is easy to collapse when the collision occurs. As for /x , /xo takes 0.14 when S is 2.15 x 10~ 2 R 
while fj,o takes 0.18 when 5 is 3.0 x 10~ 2 i?. It can be seen that roughness of the surface make 
the value of /Jo increase. 

Random lattice model can reproduce experimental tendency in (3 Q and /x with roughness 
on the surface. However, the random lattice model cannot reproduce the tendency that (3 
decreases from the maximum value f3$ in the large cot 7 ^| . Other mechanisms like sticking 
or plastic deformation on the surface may be important in the large 7 region. 

From FigEl the triangular lattice disk seems to be inadequate to reproduce experimental 
tendency. In triangular lattice, the shape of the disk is much more like a polygon than a 
circle. It may be that polygonal property of the structure causes the slip motion of the disk. 

The decrease of e in the small cot 7 in FigElcan be understood as follows. In our situation, 
normal component of initial velocity is fixed to 0.1c. Thus, the initial kinetic energy of the 
disk becomes larger in the small cot 7. As a result, the surface of the wall cannot seize the 
disk so that the initial kinetic energy is easy to propagate in the horizontal direction to the 
surface of the wall. 

In the last section we investigated the dependency on Poisson's ratio with the aid of the 
square lattice model. In contrast, we change the value of k a of the triangular lattice disk 
from 1.0 x mc 2 /R 2 to 1.0 x 10 2 mc 2 / R 2 to investigate the dependency on Young's modulus. 
Although Young's modulus increases by 100-fold, the triangular lattice disk remains slippery 
on the surface. From this fact, it can be seen that Young's modulus as well as Poisson's 
ratio are not crucial. In addition, we change the value of kf, of triangular lattice disk from 
1.0 x 10 -3 mc 2 /.R 4 to 1.0 x 10 _1 mc 2 / R 1 and investigate the effect of the nonlinear term of 
eq.(jlj). The change of kb also does not affect the results of triangular lattice disk. It can 
be seen that the nonlinear term of eq. (jH) only strengthen the surface of the model and does 
not make the triangular lattice disk rotate after collision. In the triangular lattice disk, the 
polygonal property of the surface of triangular lattice may affect the results. 

Finally, we refer to the connection between /i in eq.© and /1 in eq.Q- In FigJH the 
value of \i is same as /xq estimated from the slope in the range < cot 7 < 2 in Fig® 
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Comparing eq.fjHJ) with eq.Q, we can derive the relation between /i and ji as 

(3 X ( mR 2 \ , . 

In the two dimensional binary collision of disks, 1 + mR 2 j I can be calculated as 3 explicitly. 
Meanwhile, (3 X /[3 Z in our system can be calculated as 3.02. Thus, /i and /io are in our case 
almost identical. 



VI. CONCLUSION REMARKS 



In this paper, we demonstrate the 2-dimensional simulation of the oblique impact. Our 
random lattice model produces the same tendency as experimental data qualitatively while 
triangular lattice model can not produce the positive value of (3q. For normal COR, e depends 
on the initial angle of incidence and decreases in the large 7 when the normal component 
of initial velocity is fixed. For (3, we compare our results with Maw's theory of the oblique 
impact. Our result is consistent with their theory especially in the large and small region of 
cot 7. 

We appreciate Y. Tanaka, Y. Hayakawa and S. Nagahiro for fruitful discussion. We also 
thank M. Doi and S. Takesue for their valuable comments. Parts of numerical computation 
in this work were carried out at the Yukawa Institute Computer Facility. This work is 
partially supported by Hosokawa Powder Technology Foundation and Inamori Foundation. 
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Fig.l : The schematic figure of a collision of a disk with a wall. 

Fig. 2 : The elastic disk and wall consisted of random lattice system. 

Fig.3 : Interaction between surface particles of the disk and the wall. 

Fig.4 : The bands of random lattice (a) before and (b) after stretch. 

Fig. 5 : The schematic figures of (a) triangular lattice disk and (b) square lattice disk. 

Fig.6 : The relation between cotangent of angle of incidence 7 and f3. Cross points 
are the results of the random lattice disk. Stars are the result when the random disk 
has no roughness. Plus points are the results of the triangular lattice disk. Dashed 
and dot-dash lines are eq.Q. 

Fig. 7 : The relation between cotangent of angle of incidence 7 and COR e. 

Fig.8 : The relation between cot 7 and (3. Cross points are the numerical results of 
the random lattice model. Solid line is the theoretical curve. 

Fig. 9 : The schematic figure of the disk and the wall. A cross in a circle represents a 
center of mass of each body. 
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APPENDIX A: THE THEORETICAL DESCRIPTION OF THE OBLIQUE IM- 
PACT 

Here, we review and rewrite a theory of the oblique impact jl, 36, 3?| for our investigation. 
Let a disk with the radius R and a rectangle with the hight 2R and the width 8R be in 
contact each other as depicted in Fig.@. They have masses M and M and their radii of 
gyration k r = Rj \[2 and k' r = R a/17/3 around their centers of mass, respectively. The 
prime denotes parameters for the rectangle. The position of the contact point C is denoted 
as (r x , r z ) = (0, —R) or (r' x , r' z ) = (0, R) which are measured from the centers of mass of 
each colliding body. Here, Vi and it, (i = x, z) are relative velocity and relative displacement 
at the contact point, respectively. We assume that both normal and tangential elements of 
the compliance are proportional to the compression. We also introduce the normal stiffness 
during compression k, and the tangential stiffness during compression hi/rf for the disk. The 
equation of motion of the displacements, thus, becomes 




-m k 



f3 x v~ 2 
(3 Z 




(Al) 



where u x = d 2 u x /dt 2 , and 



, 2 



/ 2 



1 + 



mr 



mr, 



Mkl 



.2 ; 



1 + 



mr 



M'k' r Mk 2 r M 'k' r 

where 1/m— 1/M + 1/M'. In our situation, f3 x and (3 Z can be calculated as (3 X ~ 3.02 and 
& = 1. 

Equation of motion (jAl|) has two characteristic frequencies, fl and u, which are expressed 

as 

lR k I R k 1 IR 

(A3) 



+ 



mr. 



,2 " 



(A2) 



tt = 



f3 z K 



UJ = 



f3 x K 



m 



rfm 




in the normal and tangential direction to the wall, respectively. One can also express them 
using t c , which is the moment when the normal velocity of compression becomes (u z (t c ) = 
0), as = n/2t c and u = (ir/2r]t c ) \/(3 x /(3 z . According to Johnson p, stiffnesses in the 
normal and tangential direction, k and n/f] 2 , can be expressed using Young's modulus E, 
the radius of punch a, and Poisson's ratio v as 

Ea , „ 2Ea 



1 - v 2 



(2 -i/) (l + i/) 



(A4) 
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Thus, 77 can be expressed only by the Poisson's ratio, 



V 



2(1-1/)" 



2-i/ 



(A5) 



Here we define the coefficient of restitution for subsequent discussion. We assume collision 
starts at t — and compression and restitution periods terminate at t — t c and t = tf, 
respectively. The coefficient of restitution is defined as 



where p z {t) is the normal impulse, v z (t) is the normal velocity, and vo is the initial normal 
velocity. This leads to p z {tf) = (1 + e*)p z (t c ) and tf = (1 + e*)t c . In the later discussion, we 
assume that the effect of the coefficient of restitution e* is obtained by changing the stiffness 
of the normal compliant element from k to n/el at t = t c . Thus, from eq. ()A3)l . the normal 
frequency increases from Q to Q/e* at t = t c . 

1. Normal components of velocity and Force 

Now, we solve eq. (|A1)1 and obtain normal and tangential components of velocity and 
force. Assuming v z (t) + u z {z) = during contact, we obtain 



Assuming that the normal frequencies during restitution is fi/e* and the initial conditions 
v z (t c ) = and v z (tf) = — e*, we reach the exact form of solution during restitution pe- 
riod. (J A 7)1 is continuous at t = t c . By differentiating these expressions, we also obtain the 
displacement u z {t) , the force F z {t) , and the impulse p z (t) as described in Table |TJ 

2. Tangential Components of Velocity and Force 

We assume that a disk sticks or slips on the surface of a rectangle during collision and 
starts sticking at t = t<i. At t = £2, the relation between F x and F z becomes \F X \ < [iF z . We 
calculate the expressions for tangential components for sticking and slipping separately. 



Pzjtf) -Pz(t c ) 
Pz(t c ) 



Vz{tf) 



(A6) 




(A7) 
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While a disk slips on the surface of a rectangle, the relation between F x and F z is \ F X \/F Z = 
\x and the tangential velocity is changed by a impulse arising from contact force. Thus, The 
change of velocity can be described as 

1 / 

where sgn(x) = +1 for x > and sgn(x) = — 1 for x < 0. 

When a disk starts sticking at t = t 2 , the tangential oscillation starts with frequency to. 
We assume v x (t) +u x (t) = 0. By solving the equation of u x (t), we can obtain the tangential 
components of displacement, velocity, and contact force as 

u x {t2) cosu;(t — t 2 ) — ^~ 1 t> x (t 2 ) sinco>(t — t 2 ) 

uu x (t2) sinu;(t — t 2 ) + v x (t-z) cosuj(t — t 2 ) (A9) 
m(3 x l oj 2 u x (t 2 ) costo(t — t 2 ) — m(3 x l ujv x {t 2 ) sinu;(t — t 2 ) t > t 2 . 



dv x /dp 3 
dv z /dp 3 



m 



f3 xV ~ 2 
& 



u x (i) = 
v x (t) = 
FJt) = 



3. Obtaining the transition time, t\ 

We think the situation that stick begins at the initial instant contact, i.e. t 2 = 0, and 
slip begins at t — t\. From ()A9j) with the condition u x (0) = and F z (t) described in the 
Table HI t x can be obtained solving the equation 

{1 v x (0) Q sincuti 

-^^—nT = 1 0<tt<t c 

iMi^ EE^ = 1 * < . < tt { } 

^^(0)^sin(^ + f(l-e^)) ^-^KV- 

Solving these equation numerically, we obtain t\. It should be noted that there are two 
conditions if t\ is greater or smaller than t c . 

The process of initial stick takes place if t\ > 0, i.e. if in the limit as t\ — > the force 
ratio between the tangential and normal component is smaller than \i. This requires 

v x (0) ^ 2 



< (All) 
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4. Three Regimes of The Angle of Incidence 



Here, we divide all region of the angle of incidence into three regimes and calculate the 
tangential component of terminal velocity of collision for each regime. 

(i) small angle of incidence: v x (0)/v z (0) < fir] 2 

In this regime, initial stick continues until t — t\ and slip terminates at t — t/. At time 
ti, the tangential component of the relative velocity is v x {t\) = v x (0) coscuti. From eq. (|A8|) . 
the terminal tangential velocity can be expressed as 



v x (t f ) = v x (0) coswti - — \p z (t f ) - Pziti)}, 



(A12) 



where p z (tf) — p z (ti) can be expressed as 



Pz(t f ) - Pz(h) = -m^e*{l + cos(^ + -(1 - e; 1 ))} 

Pz 4 



from Table |TJ Dividing eq. (jA12|) by — v x (0) leads to 



(3 = — cosu;ti(7) — fJ-^r-e 

Pz 



1 + COS 



e 2 l 



cot 7, 



(A13) 



(A14) 



where (3 is — v x (tf)/v x (0) and cot 7 is v z (0)/v x (0). 



(ii) intermediate angle of incidence: fir] 2 < v x (0)/v z (0) < fi(l + e*)P x /P z 

In this regime, the disk initially slips and begins to stick at t = t%. After the period 

of sticking, the disk begins to slip again at t = t 3 . In the period t < t 2 , the tangential 

component of relative velocity is written as 



v x {t) = v x (0) 



m 



■Pz(t). 



(A15) 



Here let us calculate t 2 and £3. At t = t 2 , subsequent sliding and stick give the same rate 
of change for the tangential force: 



lim 



dF x (t 2 + e) 



lim ti 



dF z (t 2 + e) 
de 



(A16) 



This is the condition which determines t 2 . To simplify this condition, we need to obtain the 
exact forms of dF x (t)/dt and dF z (t)/dt. 

For tangential components of force, if one differentiate eq. ()A9j) by t, we obtain 



dF x (t) mu 3 u x (t 2 ) . mu 2 v x (t 2 ) 

- smu(t — t 2 ) cos uj{t — t 2 ) 



dt 



X 



t > t 2 . 



(A17) 



17 



Here, v x (t) is represented by eq. ()A15|) and with the aid of Table HJ and u x {t) is obtained 
from v x (t) + u x (t) = 0. Thus, the explicit expressions are 



v x (h 



u x (t 2 



v x {0) -/x^ 2 (0)[l-cosnt 2 ] 



v x (0)-fi%v z (p) 



/3 ^#sinfit 2 



{3 Z LU 



1 _ cos [ + 1(1 _ e -i; 



t 2 < t c 



to < t r 



t 2 > t c 



(A18) 



// ^o) sin fe + f(1 _ e -i ) ) , 2>tc 



For normal components of force, by differentiating the expressions of normal components 
in Table HJ we can obtain 



dF z (t 2 
dt 



-(3 z l Vt 2 mv z ($) cosfit 2 



h < t r 



n 2 mv z (0) 



COS 



fte, e* 1 2 

From (lAT7|l . (lAT8ll . and (jA"l9jl . ea. (lA"ToT) leads to 



f^ + Kl-e; 1 )) t 2 >t c 



(A19) 



fit. 



arccos 



vMhwM - p x /p 2 
v 2 - 



1>*(0) - V 



(A20) 



■ — (1 — e* ) + arccos — , — - — 



vx{0) (3x 



In the period t 2 < t < t 3 , the velocity and the force are expressed as eq. (|A9J) . This period 
of stick terminates and slip begins at time t = t 3 . t 3 can be determined by the condition 
\F X \/F Z = \i. From eq. ([A9j) and Table HI this condition leads to 



tou x {t 2 ) ttv x {t 2 ) . , , 
cos uj{t 3 — t 2 ) 77^- sin uj{t 3 — t 2 ) 



rj 2 sin 



Qt 3 7T . 1 , 



(A21) 



/i^(0) " " " w/iu*(0) 
Solving this equation numerically, we obtain t 3 . The final tangential velocity is expressed as 



v x (tf) =v x (t 3 ) - n(3 x m 1 [pz(t f ) -p z {h)], 
where p z (t) is expressed in Table [fl Dividing (jA22j) by — v x (0), we obtain 

(3 = — cosco>(t 3 — t 2 ) — [i^-[cosuj{t 3 — t 2 ) — cosfit 2 cosu{t 3 — t 2 ) 

Pz 

H sin VLt 2 sin tu(t 3 — t 2 ) + e + cos Qt 3 ] cot 7. 



(A22) 



(A23) 



(iii) Large angle of incidence: v x (0)/v z (0) > + e*)f3 x / f3 z 
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In this regime, slip does not cease before separation: ti > tf. At separation, the tangential 
velocity v x (tf) if as follows: 

v x (t f ) = v x (0) + nP x mr x (\ + e*)p g (t c ), (A24) 

where p z (t c ) = —mv z (0)//3 z . Thus, dividing ()A24|) by — v x (0), we obtain 

P = -1+^(1 + e) cot 7. (A25) 

APPENDIX B: POISSON'S RATIO OF THE SQUARE LATTICE SYSTEM 

In this Appendix, we derive the relation between elastic constants in continuum limit and 
spring constants of a two-dimensional square lattice with nearest neighbor coupling k\ and 
next nearest neighbor coupling k 2 - The elastic tensor C^ki for the two dimensional square 
lattice is represented as 

C X xxx ^yyyy ^1 ^2? (^-Q 

C XX yy Cyyxx C X yyx Cyxxy C X yxy Cy X yx ^2 5 (^^0 

and the other coefficients are zero [4^. 

Using the elastic tensor C^m and the strain tensor Uij and u^i, The free energy of the 
system U is represented as 

u = ^CijkiUijUM. (B3) 



Thus, we obtain the stress tensor 



dU 



&ij = = Cij k iu k i. (B4) 

Now we introduce the unit vector n in the axial direction of the rod. When we pull both 
sides of the rod with the pressure p, the relation 

Oik = priiUk (B5) 

holds. 

From Eqs (|B1|) . (jB2|) and (jB5|) . the explicit expressions of the stress tensors become 

&XX CxxXX^XX C xx yyUyy [k\ ~\~ k^Uxx ~\~ ^2^J/y P^X (^^) 

®yy C yyxx^xx ~t~ ^ yyyy^yy k c 2%l X x ~t~ (^1 ~t~ ky^Uyy P^y (-^^) 

@ xy ^xyxy^xy ~t~ ^ xyyx^yx 2/u2'U- C y pTL x Tly. (B8) 
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From these equations, we obtain the expressions of the strain tensors, 



u = + W ~ (B9 ) 
xx P h(k 1 + 2k 2 ) [m) 



(fa + k 2 )n 2 y - k 2 n\ 
k 1 (k 1 + 2k 2 ) 



"yy = P ~ u,u\u, ~ ( B1 °) 



u xx = p^. (Bll) 
The strain in the direction of n is expressed as u = UikUiUk- Thus, we have 

u = u xx n 2 x + 2u xy n x n y + u yy n 2 y (B12) 

h + k 2 k x - 2k 2 2 2 \ 

+ u n * n v \ V- (B13) 



k 1 (k 1 + 2k 2 ) hk 2 1 
Therefore we obtain Young's modulus as 



1 ki + k 2 k x - 2k 2 2 2 
E h(h + 2k 2 ) ^ hk 2 x y 



On the other hand, the Poisson's ratio v is defined as the ratio of the normal strain to 
the vertical strain. The latter is described as 



u t = u xx n 2 y - 2u xy n x n y + u yy n 2 x (B15) 

p. (B16) 



2(fci + 2k 2 )n 2 x n 2 y - k 2 n 2 x n 2 y 
k 1 (k 1 + 2k 2 ) k 2 

Thus, Poisson's ratio is given by 



k\ + (k 2 - ggrgrj 
A; 2 (A; 1 + A;2) + (A; 1 2 -4A;2 2 K n 2- ^ J 
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FIG. 1: H. Kuninaka and H. Hayakawa 
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FIG. 2: H. Kuninaka and H. Hayakawa 
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FIG. 3: H. Kuninaka and H. Hayakawa 
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FIG. 4: H. Kuninaka and H. Hayakawa 
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(b) 

FIG. 5: H. Kuninaka and H. Hayakawa 
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FIG. 6: H. Kuninaka and H. Hayakawa 
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FIG. 7: H. Kuninaka and H. Hayakawa 
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FIG. 8: H. Kuninaka and H. Hayakawa 
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FIG. 9: H.Kuninaka and H. Hayakawa 
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TABLE I: Normal displacement, velocity, force, and impulse 



Quantity Compression (0 < t < t c ) Restitution (t c <t <tf) 

Displacement u z (t) = -fi -1 t; z (0) sinfit u z (t) = -elVL~ l v z {ti) sin (§7 + f (l — e" 1 )) 

Velocity v x {t) = v z {G)caaQt v z {t) = elv z (0) cos + § (l - e" 1 )) 

Force F z (t) = -!^M s in^ > F z (t) = sin + |(i _ e -i)) > 

Impulse = -gg (1 - costa) p 2 (t) = -gg (l - e» cos (g + f (l - e- 1 ))) 
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